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We model sympathetic cooling of ground-state CaF molecules by ultracold Li and Rb atoms. The 
molecules are moving in a microwave trap, while the atoms are trapped magnetically. We calculate 
the differential elastic cross sections for CaF-Li and CaF-Rb collisions, using model Lennard-Jones 
potentials adjusted to give typical values for the s-wave scattering length. Together with trajectory 
calculations, these differential cross sections are used to simulate the cooling of the molecules, the 
heating of the atoms, and the loss of atoms from the trap. We show that a hard-sphere collision 
model based on an energy-dependent momentum transport cross section accurately predicts the 
molecule cooling rate but underestimates the rates of atom heating and loss. Our simulations 
suggest that Rb is a more effective coolant than Li for ground-state molecules, and that the cooling 
dynamics are less sensitive to the exact value of the s-wave scattering length when Rb is used. Lfsing 
realistic experimental parameters, we find that molecules can be sympathetically cooled to 100/rK 
in about 10 s. By applying evaporative cooling to the atoms, the cooling rate can be increased and 
the final temperature of the molecules can be reduced to 1 fjK within 30 s. 

PACS numbers: 37.10.Mn, 34.50.Cx 


I. INTRODUCTION 

Ultracold molecules are important for several applica¬ 
tions in physics and chemistry. Cold molecules have al¬ 
ready been used to test theories that extend the Standard 
Model of particle physics, for example by measuring the 
electron’s electric dipole moment mm or searching for 
changes in the fundamental constants mm- The preci¬ 
sion of those measurements can be improved by cooling 
the molecules to far lower temperatures mm- A lat¬ 
tice of ultracold polar molecules makes a well-controlled 
many-body quantum system where each particle inter¬ 
acts with all others through the long-range dipole-dipole 
interaction. This array can be used as a model system 
to study other strongly-interacting many-body quantum 
systems whose complexity is far too high to simulate on 
a computer [7]. Ultracold polar molecules offer several 
advantages for storing and processing quantum informa¬ 
tion 011], notably strong coupling to microwave photons 
and, through dipole-dipole interactions, to one another. 
The availability of ultracold molecules will also open up 
opportunities for studying and controlling chemical reac¬ 
tion dynamics in a whole new regime m- 

Some species of ultracold polar molecules can be pro¬ 
duced by association of ultracold atoms, either by pho¬ 
toassociation [II1II2] or by magnetoassociation through a 
Feshbach resonance [SHin]. Often though, the molecules 
of interest cannot be formed this way, and then more di¬ 
rect cooling methods are needed. Molecules have been 
magnetically trapped at temperatures of about 0.5 K 
by buffer-gas cooling with cryogenic helium ns HU- 
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Molecules in supersonic beams have been decelerated to 
rest and then trapped electrically and magnetically, typ¬ 
ically with temperatures in the range 1-50 mK [Mils]- 
Recently, laser cooling has been applied to SrF [3S [H] , 
YO [22] and CaF [23], and a magneto-optical trap of 
SrF has been demonstrated, producing molecules at a 
temperature of a few mK and a density of 4000 cm“^ 
[211125] • It is likely that higher densities will be reached 
using more efficient loading methods, and lower tempera¬ 
tures may be reached if sub-Doppler cooling mechanisms 
are effective. 

A promising method to cool molecules to lower tem¬ 
peratures is sympathetic cooling through collisions with 
ultracold atoms. The main difficulty with this approach 
is that static electric and magnetic traps can only confine 
molecules in weak-field seeking states, but the lowest- 
energy state is always high-field-seeking. It follows that 
inelastic collisions can heat the molecules or can trans¬ 
fer them from trapped to untrapped states. This obser¬ 
vation has motivated experimental |2hj and theoretical 
work [26I132j to search for atom-molecule systems where 
the ratio of elastic to inelastic collision cross sections is 
large. However, for most systems of interest, this ratio is 
too small for sympathetic cooling to work well. Notable 
exceptions are the Mg -|- NH system [28[, and the use of 
ultracold hydrogen as a coolant ^2], but the experimen¬ 
tal realization of those systems is exceptionally challeng¬ 
ing. An alternative approach is to use a dynamic trap, 
which could be an alternating current (ac) trap, an op¬ 
tical dipole trap, or a microwave trap, so that molecules 
can be confined in their lowest energy states. In this case, 
inelastic collisions can only excite the molecule, but the 
energy available in the collision is typically too small for 
that and so all inelastic channels are energetically inac- 
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cessible. 

In previous work [33] sympathetic cooling of a cloud of 
LiH molecules by ultracold Li atoms was simulated us¬ 
ing a very simple model. The scattering was assumed to 
be isotropic, corresponding to either s-wave scattering or 
classical collisions of hard spheres. This is appropriate 
for collisions at very low energy. However, the differen¬ 
tial cross sections at higher collision energies are typically 
peaked at low deflection angles, because many collisions 
sample mainly the long-range attraction. In the present 
work, we introduce a new collision model that takes ac¬ 
count of the full energy dependence of the differential 
cross sections. We show that this model produces sig¬ 
nificantly slower sympathetic cooling in the early stages 
than the original energy-independent hard-sphere model. 
We also consider approximations to the full model and 
show that a model that uses hard-sphere scattering based 
on the energy-dependent transport cross section [31] 
produces accurate results for the cooling of the molecules 
but not for heating and loss of the coolant atoms. 

The previous modeling work [33] explored sympathetic 
cooling in three different types of trap: a static electric 
trap, an alternating current (ac) trap, and a microwave 
trap. A static electric trap can confine molecules only 
in rotationally excited states, and it was found that for 
Li-|-LiH the ratio of elastic to rotationally inelastic colli¬ 
sions was too small for such molecules to be cooled before 
they were ejected from the trap. An ac trap can confine 
molecules in the rotational ground state, so there are 
no inelastic collisions, but elastic collisions can transfer 
molecules from stable to unstable trajectories and it was 
found that this eventually causes all the molecules to be 
lost. A microwave trap [35[ [36] can confine molecules 
in the absolute ground state, around the antinodes of 
a standing-wave microwave field, and sympathetic cool¬ 
ing in this trap was found to be feasible on a timescale 
of 10 s [33]. The microwave trap brings the benefits of 
a high trap depth and large trapping volume for polar 
molecules, especially compared to an optical dipole trap. 
In the present work, we simulate sympathetic cooling 
in a microwave trap in detail. We consider the follow¬ 
ing specific, experimentally realistic, scenario. Cold CaF 
molecules are produced either in a magneto-optical trap 
[24l [25] or by Stark deceleration [ST] [38] . In the first case 
the temperature might be about 2 mK, and in the second 
about 30 mK. The molecules are loaded into a magnetic 
trap, and then transported into a microwave trap. Here, 
the molecule cloud is compressed in order to improve the 
overlap with the atomic coolant, and this raises the ini¬ 
tial temperature of the molecules to 20 mK and 70 mK 
respectively. A distribution of atoms, either ^Li or ®^Rb, 
with an initial temperature of 100 /xK, is trapped mag¬ 
netically and is overlapped with the cloud of molecules. 
We simulate the way in which elastic collisions reduce the 
molecular temperature towards the atomic temperature. 
Black-body heating out of the rovibrational ground state 
can be reduced below 10“^ s“^ by cooling the microwave 
trap to 77K [33] . 


We start by describing our scattering calculations and 
the cross sections we obtain. Then we describe the sim¬ 
ulation method we use, and study how the choice of col¬ 
lision model affects the simulation results. Next, we ex¬ 
amine the cooling dynamics and evaluate which coolant, 
Rb or Li, is likely to be the best in practical situations. 
Because the cross section is very sensitive to the exact 
form of the atom-molecule interaction potential, espe¬ 
cially at low energies, we study sympathetic cooling for a 
range of typical values of the s-wave scattering length. In 
addition to cooling the molecules, collisions either heat 
the atoms, raising the final temperature, or eject atoms 
from the trap, reducing the atomic density. These ef¬ 
fects are particularly important if the atom number does 
not greatly exceed the molecule number. We study these 
effects and explain the results in terms of appropriate 
partial integrals over differential cross sections. Finally, 
we investigate how evaporative cooling of the atoms can 
be used to speed up the sympathetic cooling rate and 
lower the final temperature obtained. 


II. SCATTERING CALCULATIONS 

Exact scattering calculations on systems as complex 
as Li-|-CaF and Rb-|-CaF are not currently feasible. 
The combination of a deep chemical well, very large 
anisotropy of the interaction potential, and small CaF ro¬ 
tational constant mean that a very large rotational basis 
set would be needed for convergence. In addition, even if 
converged results could be achieved, uncertainties in the 
potential surface mean that no single calculation could 
be taken to represent the true system and many calcula¬ 
tions on many surfaces would be needed to explore the 
range of possible behaviors [30]. Instead we model the 
interactions with a simple single-channel model poten¬ 
tial which we choose to be the Lennard-Jones potential, 
Vir) = —C^jr^ + C' 12 /?'^^: where r is the intermolecu- 
lar distance. We have shown previously m that, while 
a simple single-channel model cannot be expected to re¬ 
produce a full coupled-channel calculation, it can quan¬ 
titatively reproduce the range of behaviors shown by full 
calculations. 

We obtain Lennard-Jones parameters for Li-|-CaF from 
ab initio calculations [33]. We obtain Cg^Li-i-CaF = 
1767 AhOQ from direct fitting to the isotropic part of 
the long-range potential, where Eh is the Hartree en¬ 
ergy and ao is the Bohr radius. We set C'i2,Li-i-CaF = 
2.37 X 10^ EhoJ^ to reproduce the depth of the com¬ 
plete potential, which is 7224 cm“^. We use the depth 
of the complete potential in preference to the depth 
of the isotropic part of the potential because the very 
large anisotropy at short-range means the isotropic part 
of the potential is not representative of the interaction. 
To obtain a Cq parameter for Rb -l-CaF we first sep¬ 
arate Ce^Li+CaF into induction and dispersion contribu¬ 
tions. Induction contributions for both systems are read¬ 
ily calculated from known values of the CaF dipole mo- 
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Center-of-mass collision energy (K) 

FIG. 1. (Color online) Total elastic cross section, (Jei (solid lines), and transport cross section, (dashed lines), for positive 
(black) and negative (red/gray) signs of the scattering length, (a) CaF+®^Rb, |a| = 1.5a; (b) CaF+®^Rb, |a| = 0.5a; (c) 
CaF+^Li, |a| = 1.5o; (d) CaF+^Li, |a| = 0.5a. 


ment [JHj and the static polarizabilities of the atoms [H] . 
The dispersion contribution for Rb+CaF can then be 
calculated from the dispersion contribution for Li+CaF 
using Tang’s combining rule [45) with known homonu- 
clear diatomic dispersion coefficients [23], atomic po¬ 
larizabilities [44] and a calculated CaF polarizability of 
aCaF = 137aQ. The sum of these contributions gives 
Ce.Rb-i-CaF = 3084i?hao- We estimate, by analogy to 
calculations on methyl fluoride [46], that the well depth 
for Rb-l-CaF will be about 2.5 times shallower than for 
Li-hCaF. This sets C'i 2 ,Rb+CaF = 1.8 x 

For our purposes, the key property of a potential is 
the s-wave scattering length, a, that it produces. In the 
present work, we vary the C 12 coefficient over a small 
range (with Cq fixed) to vary the scattering length. We 
focus on four typical scattering lengths, a = —1.5d, 
—0.5a, 4-0.5a, 4-1.5d, where a is the mean scattering 
length of Gribakin and Flambaum [57], a = 20.2 A for 
Li4-CaF and 35.7 A for Rb4-CaF. 

Discussions of thermalization have usually assumed 
that the relevant cross section is the elastic cross section 
CTei, which is the unweighted integral of the differential 
cross section da/duj, 


(Tel — 


— sin0c?0, 
doj 


( 1 ) 


where duj is an element of solid angle and 0 is the deflec¬ 
tion angle in the center-of-mass frame. However, small- 
angle scattering contributes fully to (Xei but contributes 
relatively little to thermalization. The transport cross 


section that takes proper account of this is [34], 

= 27r f ^(1 — cos0)sin0d0. (2) 

' J duj 


In the present work, scattering calculations are carried 
out using the MOLSCAT package [IS]. We use the DCS 
post-processor [53] to calculate differential cross sections, 
and the SHE post-processor [SD] to calculate 

The calculated elastic and transport cross sections for 
Li-bCaF and Rb-bCaF are shown in Fig.[^ for a vari¬ 
ety of scattering lengths. At low energy, in the s-wave 
regime, the cross sections have constant limiting values 
of 47r|ap. This is the same for both ae\ and cr^^\ be¬ 
cause pure s-wave scattering is isotropic. The cross sec¬ 
tions for positive and negative scattering lengths go to 
the same low-energy limit. However, as energy increases, 
the cross sections all diverge from one another. Those for 
negative scattering lengths, especially a = —0.5a, show 
dramatic Ramsauer-Townsend minima as the scattering 
phase shift, and hence the s-wave cross section, passes 
through a zero m- For this minimum is further 
deepened by destructive interference between s-wave and 
p-wave scattering |34j . For a = 4-1.5a a peak in both 
cross sections is seen (near 10“^ K for Rb-bCaF). This 
is a d-wave feature corresponding to the energy of the 
centrifugal barrier maximum. At higher energies, there 
are various shape resonances present for all cases. Nev¬ 
ertheless, once many partial waves contribute, the cross 
sections become less dependent on scattering length and 
approach classical limits. 
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It may be noted that the cross sections for the two 
systems for the same value of a/a are very similar, apart 
from constant factors in energy and cross section. In fact 
they would be nearly identical if the cross sections were 
in units of and energy in units oi E = 

[34l [52], where E = 9.51 mK for Li+CaF and 0.543 
mK for Rb+CaF. This scaling means that, while the 
Rb+CaF cross sections are almost independent of scat¬ 
tering length at 10 mK and above, the Li-|-CaF cross 
sections are highly sensitive to scattering length at any 
energy below 100 mK. 

For stationary atoms the molecular kinetic energy in 
the laboratory frame, E^^p, is related to the collision 
energy in the center-of-mass frame, by E^^p = 

where ^ = mcaFmat/(?7icaF + Wat) is 
the reduced mass of the collision system, mcaF is the 
molecular mass and rUat is the atom mass. The ratio 
E^^^p/E^^ is 9.40 for Li-pCaF and 1.68 for Rb-fCaF. 
This introduces a further energy scaling between the two 
systems in addition to the difference in E. 

Because the molecules are in the ground state, and the 
rotational excitation energy is far greater than the avail¬ 
able collision energy, we assume that there are no inelas¬ 
tic collisions. It is known that there can be molecule- 
molecule inelastic collisions in the presence of the mi¬ 
crowave field, even when the microwave frequency is well 
below the hrst rotational resonance [53l l54] . This is a 
concern for evaporative cooling of molecules, but less so 
for sympathetic cooling, where the density of molecules 
can be low. It is worth studying whether there can 
be atom-molecule inelastic collisions induced by the mi¬ 
crowave field, but that is beyond the scope of this paper. 


III. SIMULATION METHOD 


We assume that ground state CaF molecules are con¬ 
fined around the central antinode of a standing-wave mi¬ 
crowave field, formed at the center of an open microwave 
cavity [36j . The interaction potential of the molecules 
with the microwave field is 


U{v) 


—AU exp 


.,2 1 


COS 


2 ( 

A r 


(3) 


where AU is the trap depth and we take AU/k-Q = 
400 mK, Wx = 16.3 mm, Wy = 15.3 mm, and A = 21.3 mm 
[55] . The initial phase-space distribution of the molecules 
is assumed to be 


') — ^0,CaF 

’ ^ (27rmcaFfcBF)3/2 


X exp 


U{t) - l/(0)-fp^/(2wCaF)' 
k^T 


(4) 


where T is the initial temperature of the molecules and 
’^o,CaF is the initial density at the center of the trap, 
which is fixed such that the total number of molecules is 
iVcaF = 10^. For most simulations, we take T = 70 mK in 


order to study sympathetic cooling from a high tempera¬ 
ture. A distribution of ultracold atoms is overlapped with 
the molecules. The atoms are in a harmonic magnetic 
trap whose depth is 1 mK. We assume that the distribu¬ 
tion of atoms in phase space depends only on their en¬ 
ergy. Therefore, at all times, the atoms have a Gaussian 
spatial distribution and a thermal velocity distribution 
with temperature Tat. They have an initial temperature 
of 100^K, an initial central density of 10^^cm“^, and 
an initial number of 10®. The corresponding initial 1/e 
radius is 1.2 mm. This initial temperature and density 
can be reached by first collecting and cooling the atoms 
in a magneto-optical trap, followed by a brief period of 
sub-Doppler cooling in a molasses before loading into the 
magnetic trap. For Rb, polarization gradient cooling is 
an effective sub-Doppler cooling mechanism, while for Li 
velocity-selective coherent population trapping in a gray 
molasses can be used [551156] . Our approximation that 
the molecules are confined only by the microwave field, 
and the atoms only by the magnetic held, is a reason¬ 
able one, though our model could be extended to use the 
complete potential of both species in the combined helds. 


For each molecule, the simulation proceeds as follows. 
We solve the equation of motion in the microwave trap 
for a time step At which is much smaller than the mean 
time between collisions. Then, using the current position, 
r and velocity, v, of the molecule, we determine whether 
or not there should be a collision as follows. The ve¬ 
locity of an atom is chosen at random from a thermal 
distribution with temperature Tat. From the atomic and 
molecular velocities we calculate the collision energy in 
the center-of-mass frame, The collision probabil¬ 

ity is T = nat(r)o'(T^'^)urAt, where is the relative 
speed of the atom and molecule, riat is the ato mic d en- 
sity, and a{E^^) is either cTei or (see Sec. IV). A 
random number is generated in the interval from 0 to 1, 
and if this is less than P a collision occurs. If there is 
no collision, the velocity of the molecule is unchanged. If 
there is a collision, the velocities are transformed into the 
center-of-mass frame, a deflection angle is determined as 
described below, and the new velocities transformed back 
into the laboratory frame. If the new total energy (ki¬ 
netic energy plus trapping potential) is sufficient for the 
atom to escape from the trap, the atom, and its energy 
prior to the collision, are removed. The change in energy 
is shared among all the remaining atoms. Otherwise, the 
atom remains in the trap and the change in kinetic en¬ 
ergy is shared between all the atoms. This algorithm is 
followed for each molecule in the distribution. The den¬ 
sity and temperature of the atom cloud are updated to 
account for the atom loss and atom heating at this time 
step, and then the simulation proceeds to the next time 
step. 

With our choice of trap depth and initial atom tem¬ 
perature, there is a small evaporative cooling effect due 
to atom-atom collisions. For Rb, over the 50 s timescale 
of our simulations, 8% of the atoms are lost and the tem¬ 
perature falls to 80 /iK. Prior to Sec. IX we neglect this 
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evaporative cooling effect in our simulations because we 
wish to isolate effects that are due to atom-molecule col¬ 
lisions. Then, in Sec. |IX[ we include atom-atom collisions 
and explore the effects of evaporative cooling. 

As we will see, the molecular velocity distributions ob¬ 
tained during the cooling process are far from thermal. 
There are some molecules that never have a collision dur¬ 
ing the whole simulation and so remain at high energy 
throughout. Almost all these molecules have a kinetic 
energy greater than 10 mK, and they disproportionately 
skew the mean kinetic energy of the sample as a whole. 
Our interest is in the molecules that cool, and so we sepa¬ 
rate the kinetic energy distribution into two parts, above 
and below 10 mK. To express how well the cooling works, 
we give the fraction of molecules in the low-energy part, 
and their mean kinetic energy, both as functions of time. 

IV. COLLISION MODELS 

In previous modeling [33], atoms and molecules col¬ 
lided like hard spheres. In this model, the momenta in 
the center-of-mass frame before and after a collision, Pc 
and p^, are related by 

Pc = Pc - 2(Pc • e)e, (5) 

where e is a unit vector along the line joining the centers 
of the spheres, given by 

e = PcVl - |bp -kb, (6) 

where Pc is a unit vector in the direction of Pc and b is a 
vector that lies in a plane perpendicular to Pc and whose 
magnitude is the impact parameter divided by the sum 
of the radii of the two spheres. For each collision, b is 
chosen at random from a uniform distribution, subject 
to the constraints b • Pc = 0 and |b| < I. 

The lines labeled (i) in Fig. show how the cooling 
proceeds for CaF -|- Rb when we use the hard-sphere 
model and choose the cross section to be independent of 
energy and equal to 47rd^ = 1.59 x 10“^® m^. The cross 
section is shown in Fig. [^a), while the cold fraction and 
the mean kinetic energy of that fraction are shown in 
parts (b) and (c), both as functions of time. As explained 
in Sec. the cold fraction is defined as the fraction with 
kinetic energy below 10 mK. The cold fraction increases 
rapidly, and that fraction thermalizes quickly with the 
atoms. After just 4 s, 85% of the molecules are in the 
cold fraction and their mean energy is within 50% of the 
100 fiK temperature of the coolant atoms. 

The energy-independent hard-sphere (EIHS) model de¬ 
scribed above is reasonable at very low energy, but it has 
three deficiencies. First, it neglects the fact that the low- 
energy cross sections are actually 47ra^, where a is the 
true scattering length as opposed to the mean scattering 
length. The true scattering length can take any value be¬ 
tween —oo and -l-oo, but is generally unknown for a spe¬ 
cific system until detailed measurements are available to 




FIG. 2. (Color online) Results of various collision mod¬ 
els: (i) hard-sphere model with energy-independent cross sec¬ 
tion 47ra^; (ii) full energy-dependent differential cross sec¬ 
tion model; (iii) hard-sphere model with (iv) hard- 

sphere model with (v) hard-sphere model with 

classical approximation to The graphs show: (a) 

Cross section versus collision energy; (b) fraction of molecules 
with kinetic energy below 10 mK versus time; (c) mean kinetic 
energy of that fraction versus time. The coolant is Rb and 
a = -|-1.5a. 


determine it. Secondly, the EIHS model neglects the fact 
that real cross sections are strongly energy-dependent, 
usually showing resonance structure on a background 
that drops off sharply with increasing energy, as shown 
in Fig. I^a). Thirdly, collisions with small deflection an¬ 
gles (forwards scattering) do not contribute efficiently to 
cooling, and the EIHS model neglects the fact that dif- 
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ferential cross sections (DCS) at higher energies tend to 
be dominated by such forwards scattering, because many 
collisions encounter only the attractive long-range tail of 
the interaction potential. 

To remedy all these deficiencies, we introduce here a 
new model that we call the full DCS model. For this we 
calculate realistic integral and differential cross sections, 
as described above, for a variety of choices of the scatter¬ 
ing length a. We use the elastic cross section crei(F’®^) 
from these calculations to determine the collision prob¬ 
ability. This cross section is curve (ii) in Fig. [^a), and 
it is smaller than in the EIHS model at collision energies 
above 8mK, but larger below 8mK. We then select a de¬ 
flection angle 0 from a random distribution that repro¬ 
duces the full differential cross section, da/doj, at energy 
ifCM^ To select a deflection angle at random from this 
distribution, we form the cumulative distribution func¬ 
tion, 

5(0) = —/ ^sin(0')d0', (7) 

O'el JQ dUJ 

select a random number r between 0 and 1, and find the 
value of 0 where 5(0) = r. 

The full DCS model is our most complete one and we 
have used it for all the simulations in the following sec¬ 
tions. Its results for the choice a = -1-1.5a are shown by 
the lines labeled (ii) in Fig. It may be seen that the 
cooling proceeds more slowly than in the EIHS model. It 
takes 14 s for the cold fraction to reach 80% and for the 
energy of that fraction to be within 50% of the temper¬ 
ature of the atoms. The slower cooling is mainly due to 
the dominance of forward scattering at higher energies. 

There are three approximations to the full-DCS model 
that are worth considering because they avoid the tabu¬ 
lation of differential cross sections and cumulative distri¬ 
butions. The first of these is to use a hard-sphere collision 
model but to take the full energy-dependent elastic cross 
section from Fig.j^a). This produces the cooling behav¬ 
ior labeled (in) in Figs. ib) and (c). It may be seen 
that this model produces cooling slightly slower than the 
EIHS model, but considerably faster than the full DCS 
model. The second and more satisfactory approxima¬ 
tion is to use a hard-sphere collision model but to take 
the full energy-dependent transport cross section 
shown as line (iv) in Fig. [^a). We label this approach 
EDT-HS. It produces the cooling behavior labeled (iv) 
in Figs. I^b) and (c). It may be seen that it models the 
cooling of the molecules very accurately, because it takes 
proper account of the reduced efficiency of small-angle 
collisions for sympathetic cooling. However, as will be 
seen in Sec. |VHl[ the EDT-HS approach does not ade¬ 
quately model heating and loss of the coolant atoms. 

It is worth exploring whether a classical calculation of 
would suffice. Unlike the elastic cross section, cr^^ciass 
is finite because the factor of 1 — cos 0 suppresses the di¬ 
vergence due to forwards scattering. We have calculated 


class tbe Lennard-Jones potentials described above, 

nOO 

'^Siass = 27r/ 6[l-cos0(6)]d6, (8) 

Jo 

where b is the impact parameter and 0(5) is the clas¬ 
sical deflection function m- We find that it is very 
well approximated by the power law = 

with the dimensionless constant A = 
4.79. This cross section is labeled (v) in Fig. da). It 
agrees well with the quantum-mechanical for 

Rb-l-CaF at high energies, as we would expect when many 
partial waves contribute. Remarkably, the temperature 
and cold fraction shown for this model in Fig.j^agree very 
well with those for model (ii), even as the temperature 
approaches 100 /rK. This is an atypical result because, 
for a = -1-1.5a, cr^^ciass within a factor of about three of 

the quantum-mechanical at all energies above 3 /iK. 
For other values of a, the two cross sections can differ by 
more than a factor of three at energies below about 2E, 
which is around 1 mK for Rb-fCaF. Note that the clas¬ 
sical approximation will be less successful for a lighter 
coolant such as Li where E is far higher. 

V. APPROXIMATE COOLING RATES 

From the transport cross sections, in Fig.j^we can 
make a useful estimate of the cooling rate of molecules as 
a function of their kinetic energy. For this estimate, we 
assume stationary atoms with a uniform density riat = 
10^^ cm“^. The cooling rate is 

j plab 

^=nata(£;CM)^A£;, (9) 

where v = {2EQ^p/mciiFy^‘^ is the speed of the molecule 
and AE is the average energy transfer for a hard-sphere 
collision. AE is given explicitly as 

AE = -( -) (1 - cos0)£;J,tF- (10) 

\mCa.F + Wat/ 

Figure shows the cooling rates obtained this way, 
which although only approximate are helpful for under¬ 
standing the numerical results presented later. For colli¬ 
sions with Rb at energies above 10 mK, the cooling rate 
does not depend strongly on the s-wave scattering length. 
This is the energy regime where the a-independent clas¬ 
sical approximation to described in Sec. |iv| is 

accurate. Due to the small reduced mass in the lithium 
case, the classical limit is only reached for temperatures 
above 200 mK, and so the cooling rate depends sensi¬ 
tively on a over the whole energy range of interest. When 
a is negative there is a minimum in the cooling rates 
corresponding to the Ramsauer-Townsend minimum in 
alj^\E'^^). For Rb, at a = — 1.5d, this minimum is near 
100 fiK, which is close to the temperature of the atoms 
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FIG. 3. (Color online) Cooling rate of molecules as a func¬ 
tion of their kinetic energy, estimated from Eq. (|^, when the 
coolant is (a) Rb and (b) Li, and for various values of the s- 
wave scattering length: a = -t-1.5a (red solid line), a = -|-0.5a 
(blue dash-dot line), a = —0.5a (green dotted line), and 
a = —1.5a (black dashed line). 


in our simulations and so will not have a significant im¬ 
pact on the thermalization. For Li, the minimum occurs 
for kinetic energies between 1 and 10 mK, and so it has 
a strong effect on the thermalization. Finally, we note 
that in the ultracold limit the cooling rate is almost an 
order of magnitude higher for Rb than for Li, reflecting 
the larger value of a for Rb -I- CaF relative to Li -|- CaF. 


VI. COOLING DYNAMICS 

Figurej^shows the evolution of the {x, Vj;) phase-space 
distribution of CaF when Rb atoms are used as the 
coolant, for the case where the s-wave scattering length 
is a = -1-1.5d. At t = 0 (Fig. |^a)), the molecules fill 
the phase-space acceptance of the trap. At later times, 
more and more molecules congregate at the trap center 
as they are cooled by collisions with the atoms. After 20 s 
(Fig.|^d)), the distribution has separated into two parts. 
The majority are cooled to the center, but there are some 
that remain uncooled. These are molecules that have 
large angular momentum around the trap center and so 
are unable to reach the center where the atomic density 
is high. At a: = 3 mm for example, the atomic density, 
and hence the collision rate, is a factor of 1000 smaller 
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FIG. 4. (Color online) Time evolution of the phase-space 
distribution of molecules in the x direction. The cooling times 
are (a) Os, (b) 2s, (c) 10s, and (d) 20s. The coolant is Rb 
and a = -|-1.5a. 
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FIG. 5. (Color online) Kinetic energy distributions after 2 s, 
10 s, and 20 s. The coolant is Rb. Left panels have a = -|-1.5a 
while right panels have a = —1.5a. 


than at the center, and so molecules at this distance are 
unlikely to collide with atoms on the 20 s timescale shown 
in the figure. These molecules can be cooled by expand¬ 
ing the size of the atom cloud, but only at the expense 
of the overall cooling rate [33] . 

Figure [^a) shows histograms of the kinetic energy dis¬ 
tribution of the molecules at three different times, 2, 10 
and 20 s, when the coolant is Rb and a = -|-1.5d. These 
are the same times as chosen for the phase space distribu¬ 
tions in Fig.]^ and the results come from the same sim- 










































ulation. The initial distribution is a Maxwell-Boltzmann 
distribution with a temperature of 70 mK, truncated at 
the trap depth of 400 mK. The distribution rapidly sep¬ 
arates into two parts, those that cool and those that do 
not. The latter are the molecules that never reach the 
trap center because of their large angular momentum, as 
discussed above. A significant fraction of molecules are 
cooled below ImK after just 2 s. After 10 s the major¬ 
ity are in this group, and after 20 s this cold fraction is 
almost fully thermalized with the atoms. We return to 
part (b) of Fig.in the next section. 


VII. SENSITIVITY TO THE SCATTERING 
LENGTH AND THE CHOICE OF COOLANT 

At low energies, cross sections are very sensitive to the 
exact form of the scattering potential, as shown in Fig. 

and cannot be calculated accurately without indepen¬ 
dent knowledge of the scattering length. In our model 
Lennard-Jones potential, the full energy-dependence of 
the cross section is determined once the s-wave scatter¬ 
ing length, a, is fixed. Here, we study how the simulation 
results change as we vary the value of a. The choice of 
coolant is also a crucial consideration, and so we compare 
the results for Li and Rb as coolants. 


A. Evolution of the kinetic energy distributions 


(a)t/ = +1.5a (b)a = -1.5a 
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Figure [^compares how the kinetic energy distributions 
evolve for two cases: a = -1-1.5a and a = —1.5a, with Rb 
as the coolant. At 2 s the two distributions are similar. 
The main difference is that the distribution extends to 
lower energies for a = -|-1.5a. The similarity is due to 
the similar cooling rates at the high energies, as shown 
in Fig.|^a), while the difference at low energy is due to 
the far higher cooling rate for a = -I-1.5a at energies be¬ 
low 1 mK (compare the solid red and black dashed lines 
in Fig.j^a)). Exactly the same trend is seen after 10s 
of cooling. Once again, the high-energy parts of the dis¬ 
tributions are very similar, but the distribution extends 
to lower energies for the a = -1-1.5a case. After 20 s the 
majority of the molecules have fully thermalized with the 
atoms and the two distributions are very similar to one 
another. 

Figure shows the corresponding histograms for the 
case of Li. Here, the cooling proceeds more slowly and 
so we have added a fourth pair of histograms showing 
the distributions after 40 s. There is a great contrast be¬ 
tween the positive and negative scattering lengths in this 
case. For a = -I-1.5a the distribution evolves in a very 
similar manner to the Rb case, but when a = —1.5a it 
takes a long time for the molecules to reach energies be¬ 
low 10 mK. This is the effect of the Ramsauer-Townsend 
minimum which reduces the cooling rate estimated in 
Fig .[^b) to 0.25 s ^ for kinetic energies near 20 mK. Be¬ 
cause the minimum is broad in energy, and there is a 


FIG. 6. (Color online) Kinetic energy distributions after 2 s, 
10 s, 20 s and 40 s. The coolant is Li. Left panels have a = 
-1-1.5a while right panels have a = —1.5a. 

large mass mismatch between CaF and Li, a collision can¬ 
not take a molecule directly across the minimum. The 
molecules have to be cooled through the minimum by 
multiple collisions, and that takes a long time. Once 
molecules have passed through this minimum, cooling to 
ultracold temperatures occurs on a similar timescale to 
the a = -1-1.5d case. 


B. Cold fraction and mean kinetic energy 

Figure [^a) shows the fraction of molecules with ki¬ 
netic energy less than 10 mK, as a function of time, for 
various values of a when the coolant is Rb. This fraction 
is entirely insensitive to a. This is because the cooling 
rate is independent of a for energies above 10 mK, as we 
saw in Fig.|^ After 5 s about 50% of the molecules are in 
this cold fraction, and after 20 s this exceeds 80%. Fig¬ 
ure Bb) shows the cold fraction versus time when the 
coolant is Li. We find a strong dependence on a in this 
case. When a = -1-1.5d, the increase in the cold fraction 
with time is similar to the Rb case. For this value of a 
there is a maximum in the cooling rate at a kinetic en- 
















































9 



FIG. 7. (Color online) Fraction of molecules with kinetic en¬ 
ergy below 10 mK as a function of time for (a) Rb, and (b) Li, 
for four different values of the scattering lengths: a = -|-1.5a 
(red), a = -|-0.5a (blue), a = —0.5a (green), and a = —1.5a 
(black). 
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FIG. 8. (Color online) Mean kinetic energy of the cold 
fraction as a function of time when the coolant is (a) Rb and 
(b) Li, and for various values of the s-wave scattering length: 
a = -|-1.5a (red), a = -|-0.5a (blue), a = —0.5a (green), and 
a = —1.5a (black). 


ergy of about 70mK (see Fig.j^b)), which happens to 
match the initial temperature of the molecules, and so 
the cooling to below 10 mK proceeds rapidly. The cold 
fraction reaches 50% after 4 s in this case. The increase 
in the cold fraction is slower for a = -1-0.5d, reaching 
50% after 16 s. The accumulation of cold molecules is 
exceedingly slow when a is negative. When a = —1.5a, 
the Ramsauer-Townsend minimum is at K'caF “ mK, 
and it takes a long time for the molecules to cool through 
this minimum. The cold fraction reaches 50% after 40 s 
in this case. When a = —0.5d, the Ramsauer-Townsend 
minimum is shifted to ~ t)ut the cross sec¬ 

tion at the minimum is a factor of five smaller, and so 
the cooling is even slower, taking 50 s to reach 50%. 

Figure]^ a) shows the mean kinetic energy of the cold 
fraction as a function of time for various values of a when 
Rb is used as the coolant. As for the cold fraction itself, 
this measure is almost independent of a. This may seem 
surprising, since the cooling rates estimated in Fig.j^a) 
show a strong dependence on a below a few mK. How¬ 
ever, the mean kinetic energy is strongly influenced by 
molecules with kinetic energies close to the 10 mK cut¬ 
off that defines the cold fraction, and at this energy the 


cooling rates show little dependence on a. We find a 
small difference in the cooling rates between positive and 
negative scattering lengths. For the positive a values 
the molecular temperature is within a factor of two of 
the atomic temperature after 10 s, while for the nega¬ 
tive a values this takes 14 s. Figure [^b) shows how the 
mean kinetic energy of the cold fraction evolves when Li 
is used as a coolant. In this case, the cooling depends 
sensitively on a. When a = -1-1.5o the evolution is sim¬ 
ilar to the Rb case. The cooling is much slower when 
a = +0.5a because the low-energy cross-section is nine 
times smaller. The cooling is even slower when a is nega¬ 
tive. This is because, in the energy region between 1 and 
10 mK, the Ramsauer-Townsend minimum greatly sup¬ 
presses the cooling rate relative to the positive a case, 
and because molecules with energies in this range have a 
strong influence on the mean. 

The fraction of molecules that are cooled below 10 mK 
depends on the initial temperature, Tj. Figure com¬ 
pares this fraction for T, = 20 mK and 70 mK, for the 
case where Rb is the coolant and a = -1-1.5a. These two 
initial temperatures correspond to temperatures of 2 mK 
and 30 mK prior to compression of the cloud in the mi¬ 
crowave trap. When T; = 20 mK, more than 99% of the 
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FIG. 9. (Color online) Fraction of cold molecules as a function 
of time for the initial temperatures of Ti = 20 mK (red solid 
line), and Ti = 70mK (black dashed line). The coolant is Rb 
and a = +1.5d. 


molecules are cold within 10 s. 


VIII. ATOM HEATING AND LOSS 


The energy transferred from molecules to atoms will 
either eject atoms from the trap, or will heat them up. As 
described in Sec. m we suppose that atoms are lost from 
the trap if their total energy exceeds 1 mK. This could 
be the actual depth of the trap, or an “rf knife” might be 
used to cut off the trap at this depth. Here, we investigate 
the heating and loss of atoms and the consequences for 
sympathetic cooling. We note that while the EDT-HS 
collision model correctly captures the molecule cooling 
dynamics when is used as the cross section, it does 
not model correctly the atom heating and loss. Here, we 
highlight the difference between these two approaches by 
comparing the results obtained from the EDT-HS model 
and the full DCS model. 

Figure [T0| shows how the heating and loss rates of the 
atoms change with time in the full DCS model and the 
EDT-HS model for the case of 10® molecules and 10® 
atoms. The two models show similar trends, so we first 
discuss these trends and then consider the differences be¬ 
tween the models. At early times the majority of the 
molecules have energies far above the atom trap depth 
and so most collisions cause atom loss, rather than heat¬ 
ing. The loss rate is high while the heating rate is low. 
Nevertheless, there is still some heating due to small- 
angle collisions with the molecules which transfer only a 
little energy to the atoms. The loss rate increases dur¬ 
ing the first second because the collision cross section and 
the atom-molecule overlap both increase as the molecules 
are cooled. As time goes on the loss rate falls because 
the molecules are cooler and there are fewer collisions 
with enough energy to kick atoms out of the trap. For 
the same reason the heating rate initially increases, but 




FIG. 10. (a) Atom heating rate per molecule and (b) atom 

loss rate per molecule for the EDT-HS model (dashed line) 
and the full DCS model (solid line). The coolant is Rb, 
a = +1.5a, and the molecules have an initial temperature 
of 70 mK. There are 10® molecules and 10® atoms. 


then decreases again as the molecules have less energy to 
transfer to the atoms. For most of the 20 s period, the 
full DCS model gives more atom heating and more atom 
loss than the EDT-HS model. Only at long times, once 
the atoms and molecules are almost fully thermalized, do 
the two models give the same results. 

Integrating the results of the full DCS model shown 
in Fig.[^ we find that the total temperature increase of 
the trapped atoms is 1.3 pK per molecule, while the total 
loss is 10 atoms per molecule. The energy deposited into 
the trapped atom cloud is only 1.8% of the initial energy 
of the molecular cloud. In this sense, the sympathetic 
cooling process is remarkably efficient. 


We now turn to how the atom heating and loss rates 
can be understood, and explain why the two models give 
different results. Whether an atom is heated or lost de¬ 
pends on the kinetic energy kick it receives in the colli¬ 
sion, as given by Eq. (101 if the atoms are assumed to be 
stationary. An atom at the center of the trap is lost from 
the trap if the energy transferred in the collision exceeds 
the trap depth, AE > i?trap- This occurs if the deflection 
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angle exceeds a critical angle ©crit given by 

( 11 ) 

At laboratory-frame energies below critical energy 
Acrit = (wcaF + )/4At) Atrap, no loss is possible, as¬ 

suming stationary atoms at the center of the trap. This 
energy is 2.63 mK for Li-l-CaF and 1.04 mK for Rb-l-CaF. 
All collisions below this energy and collisions above this 
energy where 0 < ©crit will not eject atoms from the 
trap, but still transfer energy and so heat the atom cloud, 
by an amount proportional to 1 — cos©. This suggests 
the possibility of defining cross sections for atom heat¬ 
ing and loss as partial integrals of the differential cross 
section. 



cos 0 


f^loss — ^TT 


/■COS0„it 

/ —acos0, 

/_ CLLU 


( 12 ) 


^heat — 




COS ©)d cos©. 


(13) 


It is convenient to write these as integrals over d cos © 
instead of sin © d© because the cos © form allows us to 
show plots in which the integrals are simply areas that 
can be estimated by eye. Note that if < Acrit then 
CTheat = because all collisions cause heating rather 
than loss. For the full DCS model these integrals must 
be evaluated numerically, but in the hard-sphere model 
the DCS are isotropic, daus/duj = crHs/(47r), and the 
integrals can be evaluated analytically to give 

CTIoss.HS “2^^ + cos ©crit)CT hS (14) 

and 

CTheat,HS = ;^(1 — cos ©crit)^CTHS- (15) 

Figure E] shows differential cross sections at two en¬ 
ergies that correspond to Ejj^p = 2 mK and 20 mK for 
Rb-l-CaF. Both full differential cross sections and those 
from the EDT-HS model are shown (solid and dashed 
black lines respectively), and the corresponding quanti¬ 
ties weighted by 1 — cos © are shown in red. The values of 
©crit at the two energies are shown as vertical lines. In¬ 
tegrals over the complete range of cos © under the black 
lines correspond to CTcI, and under red lines correspond 
to ; the latter is the same for the full DCS and EDT- 
HS models by construction, uioss is the area under the 
black lines to the left of ©crit, and CTheat is the area under 
the red lines to the right of ©crit- It can be seen that 
at 20 mK the full DCS has a very large forwards peak; 
this dominates Uheat, even though its contribution is sup¬ 
pressed by the I — cos © weighting. The resulting CTheat 
is many times larger than in the EDT-HS model, which 
has no forward peak. The full DCS also has a secondary 



FIG. 11. (Color online) Differential cross sections and their 
contributions to heating and loss. The solid black line shows 
the full quantum-mechanical da/doj, while the solid red line 
shows (1—cos 0)da/(kjj-, the dashed lines show the correspond¬ 
ing quantities for the EDT-HS model. The vertical line shows 
the value of ©crit. (a) = 2mK. (b) = 20mK. The 

coolant is Rb and a = -|-1.5a. 


peak near cos© = 0.75, which is outside ©crit and so 
contributes to atom loss; the resulting aioss is also larger 
than in the EDT-HS model. At the lower energy of 2 mK, 
©crit is near © = tt/2. There is still a large forwards peak 
but it no longer dominates due to the changed range of 
integration, leading to similar cross sections for the two 
models. 

Figurej^shows how the heating and loss cross sections 
vary over the range of energies relevant to the cooling pro¬ 
cess. As explained above, at low energy, E^^p < Ecrit, 

we have Uheat = and ctioss = 0. Above Acrit the heat¬ 
ing cross section falls off rapidly; for the EDT-HS model 
it falls to negligibly small values by a few mK. The cross 
section for the full DCS is several times larger than that 
for the EDT-HS model in this tail, but it also falls to¬ 
wards zero. The loss cross sections for the two models 
agree surprisingly well (± ^ 30%) in an intermediate en- 
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FIG. 12. (Color online) Loss (red) and heating (blue) cross 
sections as a function of CaF laboratory energy for the EDT- 
HS model (dashed lines) and the full DCS model (solid lines). 
(7ei (solid black line) and (dashed black line) are shown 
for comparison. The coolant is Rb and a = +1.5a. 


ergy range from about 2mK to 60 mK; the extent of this 
similarity is greatest for this particular scattering length 
[a = +1.5a), but it also exists up to about 20 mK for the 
other scattering lengths investigated. Above this inter¬ 
mediate range, ctioss for the full DCS model does become 
larger than for the EDT-HS model, as we expect. The 
large peak around 1.5 mK in the elastic cross section is a 
d-wave feature that causes a large amount of backwards 
scattering around that energy; this significantly enhances 
the loss cross section because at this energy ©crit is still 
near backwards scattering. 

The overall effect is that the full DCS model gives sig- 
nihcantly larger rates of both atom heating and atom 
loss than the EDT-HS model, especially at higher ener¬ 
gies, exactly as we see in Fig.[^ This is at first sight 
surprising because each atom-molecule collision causes 
either atom heating or atom loss. However, at higher 
energies the total collision rate is considerably greater in 
the full DCS model than in the EDT-HS model, because 
the former is determined by del and the latter by 

The effects of atom heating and loss will, of course, be 
most significant when the atom number does not greatly 
exceed the molecule number. Table |T] shows the results of 
simulations for a variety of molecule numbers, with the 
atom number hxed at 10®, and once again compares the 
full DCS and EDT-HS models. In the first three rows, 
the trap depth for the atoms is 1 mK. When the atom 
number is 100 times the molecule number, atom heating 
and loss are not significant effects. For each molecule, the 
first few collisions carry away most of the energy, and al¬ 
most all of these collisions cause atom loss, rather than 
heating. Thus, for this case, 11% of the atoms are lost, 
and the atom cloud heats up by just 13 p,K. The molecules 
thermalize completely with the atoms, and the majority 
are in the cold fraction. When the atom number is only 


TABLE 1. The effect of different molecule numbers (Amoi), 
with atom number fixed at 10®, for two values of the trap 
depth Strap: 1 mK and 5 mK. The columns give the fraction of 
remaining atoms /at, the atomic temperature Tat, the fraction 
of cold molecules /moi, and the molecular temperature Tmoi 
after 50 s. The main values are for the full DCS model, and 
the values in brackets are for the EDT-HS model. 


-f^trap 

Nmol 

/at 

(%) 

Tat 

(pK) 

/mol (%) 

Tmol 

(mK) 


10'^ 

89 

(92) 

113 

(107) 

89 

(89) 

113 

(108) 

ImK 

F- 

o 

X 

lO 

38 

(59) 

159 

(136) 

88 

(88) 

168 

(144) 


10® 

2.2 

(18) 

259 

(180) 

70 

(83) 

596 

(246) 


10'^ 

95 

(96) 

151 

(133) 

90 

(89) 

153 

(134) 

5 mK 

SxlO’’ 

75 

(79) 

396 

(291) 

90 

(91) 

435 

(299) 


10® 

50 

(57) 

704 

(518) 

85 

(87) 

927 

(624) 


10 times the molecule number, the effects are far more 
dramatic. At the end of the simulation only 2.2% of the 
atoms remain, and the temperature of those remaining 
has increased to 259 /iK. Since there are so few atoms 
remaining, only 70% of the molecules now reach kinetic 
energies below 10 mK, and the temperature of this frac¬ 
tion is increased to 596 ^K. The EDT-HS collision model 
underestimates the atom loss and atom heating, and it 
predicts more cold molecules, with a lower final temper¬ 
ature, than the full DCS model. 

It is interesting to explore whether the atomic trap 
depth of 1 mK used in the simulations above is optimum. 
The last three rows of Table U show the results of sim¬ 
ulations with the atomic trap depth increased to 5 mK. 
As expected, this results in less atom loss and more atom 
heating. The fraction of cold molecules increases a little, 
but the temperature of the cold fraction increases signifi¬ 
cantly. This is especially evident when the atom number 
is only 10 times the molecule number. It is clear that 
large atomic trap depths are not necessarily beneficial 
for sympathetic cooling, and indeed there might be ad¬ 
vantages in adjusting the trap depth as cooling proceeds. 


IX. THE EFFECT OF EVAPORATIVE 
COOLING 

Evaporative cooling can be used to reduce the tem¬ 
perature further. It seems most efficient to apply the 
evaporation to the atoms, and sympathetically cool the 
molecules, rather than to apply the evaporation directly 
to the molecules. Therefore, we suppose that the evap¬ 
oration is done in the magnetic trap by applying an rf 
field which induces transitions between trapped and anti¬ 
trapped Zeeman states at a value of magnetic field only 
reachable by the most energetic atoms (an “rf knife”). 
We study the sympathetic cooling of CaF when this 
evaporative cooling is applied to Rb, for the two cases 
a = -|-1.5d and a = —0.5d. As the molecules cool, 
the molecular cloud shrinks: by choosing an appropriate 
evaporative cooling ramp, the size of the atom cloud can 












13 


be optimized throughout the sympathetic cooling pro¬ 
cess. 

We follow the theory and notation of evaporative cool¬ 
ing detailed in m- For simplicity, we assume that the 
atoms are held in a harmonic trap. The rf knife is set 
so that an atom is lost if its energy exceeds tjUbT, where 
r] is set quite large so that only the high-energy tail of 
the distribution is cut off. The rate of change of atom 
number A^at follows 

dN^t N^t 

wfw =-’ 

at Tev 

where 1/rev is the evaporation rate. It is given by 


^ev — 

V 


(17) 


where Tei is the mean time between atom-atom elastic 
collisions at the trap center. This scales with atom num¬ 
ber as 


'^el,i \ at,i / 


(18) 


where a = 77/8 — 1 and the subscript i denotes the initial 
value. Using Eqs. ( [T^ , 0 and ( |l 8 | , we obtain 


^ / iVat 

iVat.i dt Tel,! V^at.i/ 


(19) 


where k = The solution to this equation is 


iVat(t) 

iVat.i 


1 — (a — 1 ) K - 




( 20 ) 


The mean time between collisions at the start of 
evaporation is Tei,i = l/{poO'^/2v) = 70.5 ms, where 
Po = 10 ^^ cm“^ is the initial density at the trap center, 
a = 8 tt X (95ao)^ is the elastic cross section of ®^Rb at 
low temperature [58], and \/2v = 0.22m/s is the mean 
relative velocity between two atoms at the initial tem¬ 
perature of 100 /iK. The temperature of the atoms scales 
as Tat/Tat.i = , while the density scales as 

nat/nat,i = (A^at/Afat.i)^“^“'^^- In our simulations, we 
change the atom number, temperature, density and ra¬ 
dius in time, according to these results. Otherwise, the 
simulation is unchanged. We stop the evaporation when 
the atoms reach 1 pK. 

Figure |13[ a) shows how the kinetic energy distribu¬ 
tion of the molecules evolves with time when ry = 5.52 
and a = -|-1.5d. At 2 s, the distribution is similar to the 
case without evaporation (see Fig.j^a)), but by 10 s there 
is a large difference. For this value of ry the atoms ini¬ 
tially cool quickly, many atoms are ejected, and the den¬ 
sity gradually increases. About half the molecules cool 
along with the atoms and these have kinetic energy below 
100 pK at 10 s. The other half remain uncooled because 
they find themselves outside the rapidly shrinking atom 
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FIG. 13. (Color online) Kinetic energy distributions at four 
different times (2 s, 10 s, 20 s and 50 s) and for two values of 
the evaporative cooling parameter: (a) p = 5.52, (b) p = 8.14. 
For comparison, a Maxwell-Boltzmann distribution at 1 pK is 
shown by a red line. The coolant is Rb and a = -|-1.5a. 


cloud. After 50 s the cold fraction is fully thermalized 
to the 1 pK temperature of the atom cloud. Figure [T^b) 
shows the corresponding evolution when ry = 8.14. In this 
case, the evaporation initially proceeds slowly, and the 
molecule distribution remains similar to the case with¬ 
out evaporation for the first 10 s. Because the atom cloud 
shrinks more slowly a larger number of molecules are cap¬ 
tured into the cold fraction, and these then cool to 1 pK 
on a 50 s timescale. 

Figure |l4](a,b) show the fraction of molecules with ki¬ 
netic energy below 1 mK, and the mean kinetic energy 
of that fraction, using a = - 1 - 1 .5d and three different 
values of ry: 5.52, 6.67, and 8.14. When ry = 8.14 the 
atom cloud cools slowly at early times, and this gives the 
molecules enough time to thermalize with the atoms be¬ 
fore the atom cloud shrinks too much. After this initial 
thermalization to the atom temperature, the molecular 
temperature follows the evaporative cooling of the atoms 
very closely. The ultracold fraction is high in this case, 
reaching 85% after 50 s. However, it takes the full 50 s 
for this fraction to reach 1 pK. For this value of ry, the 
atom density increases by a factor of 70 over 50 s, and the 
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FIG. 14. (Color online) Sympathetic cooling of molecules 
with evaporative cooling applied to the atoms. Graphs show 
the time evolution of (a) the fraction of molecules with kinetic 
energy below 1 mK, when a = +1.5a; (b) the mean kinetic en¬ 
ergy of the ultracold fraction when a = -|-1.5a; (c) the mean 
kinetic energy of the ultracold fraction when a = —0.5a. (i, 
black) 77 = 5.52, (ii, red) 77 = 6.67, (hi, blue) 77 = 8.14. In (b) 
and (c), the dashed lines show how the atomic temperature 
evolves. The long-dash green line shows the atom tempera¬ 
ture without evaporative cooling. 


mean atom-molecule collision rate increases from 4 s“^ to 
45s“^. When 77 = 6.67 the atoms cool more rapidly and 
the cloud size shrinks more rapidly. Consequently, the ul- 
tracold fraction of molecules is reduced to 74% but this 
fraction now reaches 1/rK in 30 s. When 77 = 5.52 the 
atoms initially cool quickly, but the cooling rate slows 
down as time goes on because the density does not in¬ 


crease rapidly enough to compensate for the decrease in 
atom velocity. The ultracold fraction of molecules re¬ 
duces to 59%. The mean kinetic energy of this fraction 
falls quickly, reaching 100 /iK in 3.4 s, and 10 /iK in 17 s. 
Therefore, evaporative cooling with a relatively low 77 is a 
good strategy for cooling rapidly to temperatures above 
10 /rK. However, the cooling slows down at longer times 
and it ultimately takes longer to reach 1 /rK than for the 
intermediate value of 77 . 

Finally, we consider the case where a = —0.5d. This 
is a highly unfavorable case compared to a = -|-1.5d, 
both because the elastic cross section in the ultracold 
limit is nine times smaller and because there is a deep 
Ramsauer-Townsend minimum in the cross section for 
collision energies slightly below 100 77 K, as can be seen 
in Fig.[^ We find that the fraction of molecules with 
kinetic energy below 1 mK is almost unchanged from that 
shown in Fig. (Ha). This is to be expected since, at 
energies higher than 1 mK, the cross sections for the two 
values of a are not too different. Figure T^c) shows how 
the mean kinetic energy of the ultracold fraction evolves 
when a = —0.5d. Because of the lower collision rate, the 
mean kinetic energy of the molecules lags behind that 
of the atoms, instead of the two being locked together 
as they are in the case of a = -1-1.5d. The molecules 
are slow to reach 20 /iK for all values of 77 , because they 
have to cool through the Ramsauer-Townsend minimum 
to do so. For 77 = 5.52, the atoms cool too quickly and 
the molecules have not thermalized with the atoms even 
after 50 s. For 77 = 8.14 the initial cooling rate of the 
atoms is slow enough that the molecule temperature can 
more closely follow the atom temperature, both reaching 
1 /iK in about 50 s. The cooling of the molecules is fastest 
for the intermediate value of 77 . In particular, the mean 
kinetic energy of the molecules falls rapidly as soon as it 
is below 20 /iK, and it reaches 1 /iK in 36 s. We see that, 
even for this unfavorable value of a, evaporative cooling 
of the atoms can bring the molecule temperature down 
to 1 /iK on a reasonable timescale, provided a suitable 
value of 77 is chosen. It is clear that knowledge of the 
actual atom-molecule scattering length will be needed to 
choose the optimum conditions for evaporative cooling. 


X. CONCLUSIONS 

In this paper, we have addressed the methodology for 
modeling sympathetic cooling of molecules by ultracold 
atoms, and we have studied in detail the results of sim¬ 
ulations for a prototype case where ground-state CaF 
molecules in a microwave trap are overlapped with ultra¬ 
cold Li or Rb atoms in a magnetic trap. This work leads 
to a number of conclusions which we now summarize. 

Previous work on sympathetic cooling used a hard- 
sphere model of collisions based on an elastic cross sec¬ 
tion. This is appropriate at very low energies (in the s- 
wave regime), but breaks down badly for heavy molecules 
in the millikelvin regime. We have shown that a hard- 
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sphere model based on an elastic cross section signifi¬ 
cantly over-estimates the cooling rate for collision en¬ 
ergies above the s-wave scattering regime. A hard- 
sphere collision model that uses the energy-dependent 
momentum transport cross section, gives the cor¬ 

rect molecule cooling rate, but underestimates both the 
heating of the atoms and the loss of atoms from the trap. 
We have therefore used the full differential cross section 
to model atom-molecule collisions, so that the cooling 
of the molecules and the associated heating and loss of 
atoms are all modelled accurately. 

We have studied sympathetic cooling of CaF with both 
Rb and Li over a range of typical values of the atom- 
molecule scattering length a. We find that Rb offers sig- 
nihcant advantages over Li as a coolant for ground-state 
molecules. The mean scattering length a is almost twice 
as large for Rb, and so it is likely that the true scatter¬ 
ing length will also be larger for Rb. The mean energy 
transfer is proportional to /i/(mcaF + mat) which is 0.48 
for Rb, but only 0.19 for Li. If a happens to be negative 
there can be a deep Ramsauer-Townsend minimum in 
the cross section. For Li, the minimum typically occurs 
when between 1 and 10 mK, and the molecules 

cool very slowly because their energies must pass through 
this minimum. For Rb, the minimum is shifted down an 
order of magnitude in energy, and so the molecules do 
not encounter the minimum until they have reached the 
ultracold regime. For Li, the cooling rate is very sensitive 
to the actual value of a, while for Rb the initial cooling 
rate is fairly insensitive to a because the Rb-|-CaF cross 
section conforms closely to a classical result, independent 
of a, down to temperatures near 1 mK. This brings less 
uncertainty about the likely results of sympathetic cool¬ 
ing experiments if Rb is used. These advantages of Rb 
as a coolant are likely to extend to other molecules of a 
similar or greater mass. Finally, it is experimentally eas¬ 
ier to prepare large, dense samples of ultracold Rb than 
of ultracold Li. 

It should be noted that the preference for Rb over Li 
applies only to ground-state molecules that cannot be lost 
from the trap through inelastic collisions. For molecules 


in static magnetic or electric traps, a light collision part¬ 
ner such as Li, Mg or H provides a higher centrifugal 
barrier than a heavy one such as Rb, and this may be 
important for suppressing low-energy inelastic collisions 

Hi [33 [35]. 

For molecules with an initial temperature of 70 mK, 
cooled by Rb with a temperature of 100 p,K and a peak 
density lO^^cm”^, we find that, after 10s, 75% of the 
molecules have cooled into a distribution with a tempera¬ 
ture of 200 /iK. If the initial temperature of the molecules 
is reduced to 20 mK, this fraction increases to 99% due 
to improved overlap between molecule and atom clouds. 
By arranging for the atom trap depth to be far below 
the initial molecule temperature, we can ensure that the 
majority of the energy in the molecule cloud is removed 
by atoms that are lost from the trap, instead of heating 
the atom cloud. For efficient cooling the atom number 
should exceed the molecule number by at least a factor of 
100. By applying evaporative cooling to the atoms, the 
molecules can be sympathetically cooled more rapidly, 
or they can be cooled to far lower temperatures. For 
values of the scattering length in the likely range, and 
with a suitable choice of evaporation ramp, 70% of the 
molecules can be cooled to 1 fjK within about 30 s. These 
are all encouraging results: using experimentally achiev¬ 
able atom numbers, densities and temperatures, sympa¬ 
thetic cooling to ultracold temperatures can work on a 
timescale that is short compared to achievable trap life¬ 
times. A good starting point for such experiments would 
be a mixed-species magneto-optical trap of molecules and 
atoms. 

Data underlying this article can be accessed at 
http://dx.doi.org/10.528I/zenodo.32993 and used under 
the Creative Commons CCZero licence. 
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